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ABSTRACT 


A unique numerical method has been developed to model the dynamic processes of bulk 
evaporation and condensation processes, associated with internal heat generation and natural 
convection under different gravity levels. The internal energy formulation, for the bulk liquid- 
vapor phase change problems in an encapsulated container, was employed. The equations, 
governing the conservation of mass, momentum and energy for both phases involved in phase 
change, were solved. The thermal performance of a multiphase uranium tetrafluoride fuel element 
under zero gravity, micro-gravity and normal gravity conditions has been investigated. The 
modeling yielded results including the evolution of the bulk liquid-vapor phase change process, 
the evolution of the liquid-vapor interface, the formation and development of the liquid film 
covering the side wall surface, the temperature distribution and the convection flow field in the 
fuel element The strong dependence of the thermal performance of such multiphase nuclear fuel 
cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10' 3 -g and 
1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall 
is more isothermalized at the wall surface, which can prevent the side wall from being overheated. 
As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the 
liquid film and smaller across the vapor phase. This investigation provides valuable information 
about the thermal performance of multi-phase nuclear fuel element for the potential space and 
ground applications. 
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Nomenclature 


A, C, D constants 

E internal energy [J] 

G internal heat generation rate [kW/kg] 

H enthalpy [J] 

O heat quantity [J] or volumetric thermal power density [MW/m 3 ] 

T temperature [K] 

V volume [m 3 ] 

U velocity [m/s] 

W work [J] 

X, R dimensionless coordinates 

c Vi c p specific heat at constant volume and constant pressure [J/kg/K] 

e specific internal energy [J/kg] 

h specific enthalpy [J/kg] 

/ vapor phase fraction 

L length scale or the length of the container [m] 

p pressure [bar] 

t time [s] 

x, r cylindrical coordinates [m] 

Nondimensional Characteristic Numbers 

Gr Grashof number 
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Pe 


Peclet number 


Pr Prandtl number 

Re Reynolds number 

Ra Rayleigh number 

St Stefan number 

Greek symbols 

a thermal difEusivity [m 2 /s], or thermal diffusivity ratio 

e phase-change temperature interval 

k thermal conductivity [W/K/m], or its ratio 

O order of magnitude 

0 dimensionless temperature 

p density [kg/m 3 ], or its ratio 

x dimensionless time 

Subscripts 

/ interface position, mixed phase 

/ interface 

/ liquid phase 

r characteristic scale 

s, sat saturation point 

v, vap vapor phase 

Head Bar 

nondimensional quantity 
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Chapter I. Introduction 


The thermionic energy converter is a non-mechanical gaseous electronic device for 
converting heat energy directly into electric potential through thermionic electron emission In a 
thermionic diode, electrons are emitted from a hot electrode, which is the emitter, and collected 
by a colder electrode, which is the collector, at a higher potential energy (lower electrical 
potential). Part of the heat removed from the emitter due to evaporating electrons is rejected to 
the collector by condensing electrons The remaining heat is converted into electric power in the 
load as electrons return to the emitter potential. The physics governing thermionic energy 
conversion is well understood (Hatsopoulos and Gyftopoulos, 1973), and this energy conversion 
technique is attracting wide interest in its application toward a variety of power conversion 
systems (Lee et a!., 1993, Young et al., 1993). In a review paper on thermionics, Rasor (1991) 
outlined the history, application options, and ideal performance of thermionic energy converters 
and described the basic plasma types associated with various modes of converter operation 

Various types of thermionic-converter systems, such as solar thermionic generators, 
radioisotope thermionic generators, chemical thermionic generators and nuclear reactor 
thermionic generators, have been or are currently being developed. The nuclear reactor thermionic 
generators are simply thermionic converters supplied with heat produced by nuclear fission. The 
fission heat is conducted from the reactor core to the emitter, and the rejected heat is conducted 
from the collector and carried away via a coolant A typical schematic configuration of the 
thermionic fuel element (TFE) fueled with solid nuclear fuel UO 2 is shown in Fig. 1.1 The nuclear 
fuel resides in the center of the TFE. The space between the emitter and collector is filled with an 
easily ionizable rarefied medium, typically cesium vapor, in order to enhance the electrical 
performance. A liquid metal coolant flows in the annulus channel surrounding the TFE. Recently 
a variety of investigations involving such a TFE have been carried out that consider the thermal 
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and electrical performance, system simulation and design, and application as a space power 
system, Lewis etal. (1991), McVey and Rasor (1992), Young etal. (1993). 
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Figure 1 . 1 Schematic of a typical nuclear thermionic fuel element 


The thermal performance of a thermionic fuel element is tightly coupled with its electrical 
performance. The electrical output current density is directly related to the emitter temperature, 
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that is, the current density increases significantly with emitter temperature, which can be 
illustrated by a simple example of an ideal diode thermionic converter as follows 

J = AT e 2 exp HV + * e )lkT E ] 

where J is the current density, A and k are constants, V is the output, f c is the collector work 
function, and T E is the emitter temperature. 

Consequently a nuclear fueled TFE, which can provide high temperature, is well suited for 
this application The maximum operating temperature of solid nuclear fuel U0 2 could be as high 
as 1800 K In order to further enhance the TFE output, a higher operating temperature is 
required. Such higher temperatures could be realized if liquid nuclear fuel is employed. If nuclear 
fuels such as UF 4 or UF 4 /UO 2 in a liquid or vapor state were used, they could sustain a 
temperature of 2600 K or higher. The large increases in operating temperature will raise the 
power output significantly The use of a multi-phase nuclear fuel reactor in conjunction with a 
thermionic converter provides the potential for a compact, higher power output electrical 
generator that is well suited for space applications. 

Certain problems related to the thermal performance of TFE arise with regard to the 
application of multi-phase nuclear fuel reactors. These include the phase change (PC) process 
between the liquid and vapor, the behavior of evaporation and condensation associated with 
internal heating and external cooling, the location and motion of the liquid-vapor interface, the 
temperature distribution, the influences of gravitational fields, and so on. With such a system, 
particular attention must be paid to limiting the vapor pressure that is safe for the containment 
vessel because evaporation in a closed container will raise the vapor pressure very rapidly. A 
thorough investigation of these issues will be required to fully apply this concept 
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Figure 1.2 Schematic of the multi-phase nuclear fuel system 

A cylindrical container filled with multi-phase nuclear fuel, such as UF 4 , is shown in 
Figure 1.2. The fuel experiences liquid-vapor phase change processes: evaporation and 
condensation. The internal heat provided by the nuclear fission reaction evaporates the liquid, 
while the side wall cooling tends to condense the vapor. Globally, the evaporation of liquid and 
condensation of vapor take place simultaneously. The evaporation and condensation processes 
could be driven by either pure heat conduction or by both heat conduction and convection, 
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depending mainly on the gravity field. Under a zero-gravity condition, buoyancy-driven 
convection is nonexistent. However, under micro-gravity or one-gravity fields, buoyancy-induced 
convection is present. The location of the liquid-vapor interface may vary depending on the mode 
of heat transfer, the internal heat generation intensity and spatial distribution, and the gravity field. 
Such an interface could be defined as a free and moving interface. To solve the phase change 
problem analytically or numerically, the free and moving interface is not known in advance, so 
that its position has to be determined as a part of the solution. Thus the solution to this heat 
transfer problem is quite complex. 

The purpose of this work was to develop a numerical method for the computation of the 
bulk evaporation and condensation processes involved in a specific class of liquid-vapor phase 
change problems under constant volume. The phase-change behavior and thermal performance of 
a multiphase nuclear fuel element, operated under 0-g, micro-g and normal-g, were investigated. 
Certain problems are related to the application of the multiphase fuel element which include the 
phase-change processes between liquid and vapor, i.e., the bulk evaporation and condensation, 
with internal heat generation, the location of the liquid-vapor interface, the temperature 
distribution, the influences of gravitational fields, etc. The evaporation and condensation 
processes in the multiphase system could driven by either pure heat conduction or by both 
conduction and convection, depending mainly or the gravity level. Under zero-gravity condition, 
buoyancy driven convection does not exist. However, under micro-gravity or normal gravity 
fields, buoyancy induced convection is present. The location of the liquid-vapor interface may 
vary depending on the mode of heat transfer, the convection flow pattern, the internal heat 
generation rate and its spatial distribution, as well as the gravity field. 
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Specific issues had to be investigated for the above desired conditions. The literature 
appears not to contain significant work on the phase change processes between liquid and vapor, 
the bulk evaporation and condensation involved in the multiphase system Very limited 
investigations have been conducted, i.e., Shyy (1994b), Ding (1995), Ding and Anghaie (1994b, 
1995a, b), in this area To the best of the author’s knowledge, no prior research has considered 
the general bulk condensation and evaporation processes in an enclosed container, under various 
gravity fields, with and without internal heat generation. Liquid-vapor phase change problems are 
also characterized by a large step change in the density across the phasic boundary. The singular 
jump of the density at the interface, other than the jump of temperature gradient which exists 
unavoidably for all phase change problems, brings a serious challenge to computational fluid 
dynamics and heat transfer 

Since the investigations for the solid-liquid phase change problems, such as melting and 
solidification, are numerous, some successful experiences and lessons can be drawn and applied 
into the study of the bulk evaporation and condensation. Appropriate modeling, such as 
formulation and numerical treatment have to be developed. For the solid-liquid phase change 
problems, the thermodynamic and transport properties between two phases vary moderately 
during phase change processes. The properties, however, vary significantly or dramatically during 
the phase change processes between liquid and vapor. Therefore, scaling and analysis of order of 
the magnitude have to be conducted to make simplifying assumptions that can lead to a concise 
and efficient solution procedure for the governing equations. Special treatment and closure 
relationships, such as those for vapor phase fraction, the vapor density, the saturation pressure 
and temperature, had to be applied and implemented into the computation. 
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Like all phase change processes, evaporation and condensation are also inevitably subject 
to the effect of convection induced by either buoyancy or surface tension variations The existence 
of the convection can substantially affect the flow field, the temperature field, and subsequently 
the performance of the phase change process. The effect of the convection on the evaporation and 
condensation processes are to be investigated. For the conventional solid-liquid phase change 
problems, heat is usually provided externally from the boundary to the system. For the nuclear 
system, heat is generated internally via nuclear fission reaction. The problems of phase change 
induced by internal heating are different from those induced by external heating and need to be 
investigated too. The internal heating, like the external heating, results in the temperature and 
density distribution, and consequently causes the buoyancy effect. Such natural convection 
induced by internal heating is defined as internal natural convection, to be distinguished with the 
classical natural convection caused by external heating/cooling. This internal natural convection 
has drawn the attention of a few investigations, i.e., Acharya and Goldstein (1985), Lee and 
Goldstain (1988), Ding(1995). The behavior of the internal natural convection and its influence on 
the phase change processes, in particular on the bulk evaporation and condensation processes, 
were investigated in these studies 

Previous investigations of phase-change problems are numerous but limited to a few types. 
Those activities, which can be traced back as early as Clausius, Clapeyron, Stefen, Nusselt, etc . , 
have been conducted both experimentally, analytically and are currently being attempted using 
numerical simulations and approximation methods. In their annual review papers, Eckter et al, 
(1990, 1991 and 1992), suggested that most of these current works have focused on: melting and 
solidification, pool and flow boiling, film and droplet evaporation, filmwise and dropwise 
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condensation, with applications in thermal energy storage, solidification of alloys and metals, 
casting processes, crystal growth, heat pipes, heat pumps, boilers, heat exchangers, and nuclear or 
chemical reactors. The major characteristics and classification of a group of phase change 
problems, and the current research activities and status are summarized in Table 1.1. 


Table 1 . 1 Brief summary of major characteristics of phase change problems 


Moving Boundary - Phase Change Problems 

PC Between 

Solid - Liquid 

Liquid - Vapor 

Process 

Melting and 
solidification 

flow or 
pool 
boiling 

film or 
drop 
conden- 
sation 

bulk evaporation 
and 

condensation 

Heating 

/Cooling 

externally from 
the system boundary 

externally 
from boundary 

externally or 
internally heating 

Thermo- 

physical 

properties 

constant or 
change moderately 

change significantly or drastically 

Conduction 

and/or 

convection 

conduction 

only 

conduction 

and 

convection 

conduction and 
convection 

conduction 

only 

conduction 

and 

convection 

Solution 

Approach 

analytical 

numerical 

empirical 

numerical 

empirical 

experimental 

numerical 

analytical 

numerical 

Research 

Status 

numerous investigations 

numerous 

investigations 

very few found, 

being conducted 
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Early analytic work on phase change between solid and liquid included studies by Lame 
and Clapeyron (1831) and Stefan (1891). The fundamental features of this type of problem and 
some basic analytical solution methods were discussed by Ozisik (1980), Crank (1984) and Yao 
et al. (1989). Crank (1984) presented a broad but reasonably detailed account of the mathematical 
solution, both analytical and numerical, of these problems. Due to the existence of the interface, 
difficulties were encountered in solving these problems to the extent that only a handful of very 
simple cases can be theoretically analyzed. 

Since the 1970’s, greater attention and interest have been paid to purely numerical analysis 
and simulations, which is reflected in the large number of publications covering these topics, i.e., 
Crank (1984), Voller et al. (1987, 1991), Yao et al. (1989), Brice (1986), Brown (1988), Gau 
and Viskanta (1986), Nadarajah and Narayanan (1990), Shyy et al. (1990, 1991, 1992a, b, 1994a, 
1995) and Shyy (1994a, b, 1995). In his book, Shyy (1994a) summarized a substantial portion of 
his research activities, and presented the various computational elements important for the 
prediction of complex fluid flows and interfacial transport. 

Recent years have also witnessed a growing interest in other forms of phase change processes, 
such as evaporation and condensation. Carey (1992) systematically summarized the state of the 
art in boiling and/or condensation phenomena. That book focuses on some fundamental elements 
and the basic physical mechanism of condensation and vaporization processes and provides a 
broad review of the research activities in these areas. Among the investigations on condensation 
and evaporation, those on film condensation are well developed, particularly on laminar film 
condensation, due to its simplicity. The development of the current theory and the prediction of 
heat transfer rates associated with film condensation have its origin in the analyses presented by 
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Nusselt (1916). Subsequently, many modifications of this analysis have appeared, each one 
further relaxing the restrictive assumptions of the Nusselt analysis. Considerations have been 
given to gravity-flow condensation, pure forced convection condensation in the presence of a 
body force, such as in Sparrow et al. (1967) and Rohsenow (1973). Recently some studies on 
general film condensation transients and the interfacial condensation model have been conducted 
by Filk and Tien (1989) and Gemer and Tien (1989). With the development of computational 
fluid dynamics and numerical heat transfer, more and more attempts have been made to model and 
simulate condensation phenomena, for example, Jones and Renz (1974), Cossmann, Odenthal 
and Renz (1982), and Gaultier et al. (1993). However, most of these activities, either analytical 
or numerical, have focused on flow condensation problems with open plates, channels or open 
containers where the vapor pressure, saturation temperature and latent heat are maintained 
constant. These problems, where forced convection plays a dominant role, are usually parabolic 
type problems, and consequently are much simpler than those in the present study in which the 
saturation temperature, pressure and latent heat continuously change. A few computational 
simulations of evaporation in containers have been carried out recently. Jang et al. (1990) 
modeled a heat pipe start-up from the frozen state by using one dimensional equations for 
different phases. Shyy (1994b) developed a computational method for predicting the two-phase 
transient flow and heat transfer characteristics within a constant pressure reservoir of a capillary- 
pumped-loop and investigated the phase change of liquid and vapor, under various conditions. 
Ding and Anghaie (1994a, 1995a, b), Anghaie and Ding (1995), Ding (1995) developed a 
modeling procedure for the bulk evaporation and condensation and applied it into the 
investigation of multi-phase nuclear fuel systems. 
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Considerable attention has been given to the study of natural convection of single phase 
fluid in enclosures. The temperature gradients across the system are supplied externally to induce 
the density gradients and consequently the buoyancy driven flows. A great deal of research, both 
numerical and experimental, on various types of convective heat transfer problems has been 
conducted, i.e., Ostrach (1982, 1983, 1988), Gebhart et al. (1988), de Vahl Davis (1983), 
Hortmann et al. (1990), Shyy et al. (1992b), Shyy (1994a), Bejan (1984), Carpenter et al. (1989), 
Ramaswamy et al. (1992), Bergman et al. (1981), etc.. Ostrach (1982, 1983, 1988), Gebhart et 
al. (1988) and Bejan (1984) reviewed and summarized those activities and achievements. Among 
the large number of the natural convection problems, the buoyancy-driven flow in a square cavity 
with vertical side walls heated differentially, so called a bench mark problem for natural 
convection, is a suitable vehicle for testing and validating computer codes. Accurate benchmark 
solutions have been published by De Vahl Davis (1983) and by Hortmann et al. (1990) and are 
available for comparison with the present computations. 

Natural convection induced by external heating is of importance in many applications 
Natural convection in an enclosure in which the internal heat generation is present is also of prime 
importance in certain technological applications. Examples are post-accident heat removal in 
nuclear reactors and some other geophysical or astrophysical problems. Natural convection 
induced simultaneously by external heating and internal energy sources, however, in contrast to 
the extensive investigations on the conventional external heating convection problems, has 
received very limited attention, such as Archarya and Goldstein (1985), Lee and Goldstein (1988), 
Fusegi et al. (1992), Ding (1995). The presence of internal heat generation provides an additional 
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dynamic in overall convective flow systems and offers a complicated yet challenging aspect of the 
convection under consideration. 

With no exception, convection induced by either buoyancy or surface tension variations will also 
be involved in evaporation and condensation like all phase change problems. The involvement of 
the convection can substantially affect the phase change process. Numerous investigations on the 
solid-liquid phase change problems with convection can be found in Ostrach (1983), Crank 
(1984), Brice (1986), Brown (1988), Gau and Viskanta (1986), Nadarajah and Narayanan (1990), 
Shyy etal. (1990, 1991, 1992a, b, 1994a, 1995) and Shyy (1994a, b, 1995), etc.. 
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Chapter II. Thermal Performance of Multiphase UF4 Fuel Element at Zero Gravity 

2. 1 The Internal Energy Formulation 

The first law of thermodynamics, for any closed thermodynamic system, is 

SO = dE + SW = dH-(pdV + Vdp) + SfV (2.1) 

where E is the internal energy, H is enthalpy, Q is the heat input, W is the work output, p and V 
are the pressure and volume of the system, respectively. For the phase change process involved in 
the rigid and enclosed container with constant volume, there is no mechanical work Therefore the 
first law for this closed and rigid system reduces to 

SQ = dE= dH - Vdp (2.2) 

which indicates that the heat exchange between the system and the ambient causes the internal 
energy, the enthalpy and pressure of the system to change. Therefore the released or absorbed 
latent heat associated with the phase change induces the changes of internal energy, enthalpy and 
pressure. The quantity of the heat exchange is equal to the internal energy change, but not equal 
to the enthalpy change. The conventional enthalpy formulation of the energy equation, commonly 
used for the phase change between solid and liquid, should include a transient pressure term for 
this phase change between liquid and vapor The weak form of the energy equation based on 
enthalpy formulation is 



-P- + V • (kVT) 
Dt 


(2.3) 


This formulation including the transient pressure term introduces additional complexity 
and does not seem to be the best choice for this problem. A more attractive alternative, the 
internal energy based formulation for the energy equation which does not include the pressure 
term, is proposed and developed in this work. Similar to the enthalpy formulation used for melting 
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and solidification problems, the internal energy formulation does not require to explicitly track the 
interface position. The weak form of the energy equation with internal energy as the primary 
variable can be expressed as follows: 


,- = V( ^ (2.4) 


where the internal energy e may be expressed, in the form of a summation of sensible heat and 
latent heat, for the liquid, vapor and mixed phases as: 


« = c ,,i T f = 0 

* = C v,l T sat + / Ag 0</< (2.5) 

e = c T +Ae + c (T-T ) f = 1 

v,/ sat v,v sat 

where T sat is the saturation temperature, Ae is the internal energy difference between saturated 
vapor and saturated liquid (latent heat), and the constant volume specific heat c v ’s are functions of 
temperature provided by Reynolds (1979). The vapor phase fraction /is zero in the region filled 
with liquid phase and unity in the region occupied by vapor phase. The vapor phase fraction lies 
between zero and unity when the control volume is undergoing phase change. 

Substituting Eq.(5) into the governing Eq.(4) yields 


A^-fe,r) = v( «r , vr» 

/ = 0 



0</<l 

(2.6) 

A^fe,0 = V-(«:,VD 

/= 1 



where the latent heat now appears as a source term. The form of this internal energy based 
equation is identical to that of the enthalpy based equation, except the pressure term appearing in 
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the later one. Just like the enthalpy formulation, this internal energy formulation is also a single 
region formulation, where one set of governing equations can describe both phases A major 
advantage of the internal energy formulation is that it is not necessary to split the domain into 
separate subdomains consisting of different phases A fixed grid can be employed to facilitate the 
computations. 

Both enthalpy formulation and internal energy formulation are implemented and examined 
The vapor phase fraction, /, may be defined in two different ways, which can lead to different 


update methods. 

Some additional closure relations are required for the multiphase condition in the 
container. Property variations are based on a simplified searching routine. The different phases are 
differentiated based on the local temperature or internal energy at each nodal point, within the 
frame work of the temperature based or internal energy based update method. Using temperature 
based update, the liquid or the vapor phase is identified by the local temperature, wherever the 
temperature at any point is lower than the saturated liquid temperature or higher than the 
saturated vapor temperature, that point is located in the liquid or vapor region, respectively. 
Similarly, if the internal energy based update is used, the liquid or the vapor phase is identified by 
the local internal energy. Either case allows the vapor phase fraction in a computational cell to be 
calculated Using the Clausius-Clapeyron equation, the saturation temperature as a function of the 
system pressure is calculated: 


C 

Tsat ~ 

D Psat 


( 2 . 7 ) 


where C and D are constants. 
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The saturation pressure, or the vapor pressure, can be obtained from the approximate 
equation of state for an ideal gas, and the mass conservation in which the bulk vapor density and 
temperature are used: 

Psa, = P vap = m total = m v + mi = Constant (2.8) 

r vap 

where R is the gas constant for vapor phase. Eqs.(2.7) and (2.8), have been proposed and used as 
the basis for the model developed by the authors. Ding and Anghaie (1994a, b, 1995a, b, c), 
Anghaie and Ding (1995) 

2.2 The Governing Equations 

Based on the modeling procedure developed in the previous section, the bulk evaporation 
and condensation processes with internal heat generation under zero gravity condition were 
studied. The system, at a certain steady state, filled with liquid and vapor is shown as Fig. 1.2. 
Since the system is designed to operate under 0-g of the space environment, there is no buoyant 
effect and natural convection Initially, the vapor and liquid are kept saturated at a certain 
pressure and uniform temperature The top wall and bottom wall are insulated. The side wall is 
cooled with a lower constant temperature or constant heat flux. The internally generated heat 
evaporates the liquid, while the cold side wall removes the heat from the liquid-phase and 
condenses the vapor-phase. For this system with such an orientation and operating under 0-g, 
there is no bubble or droplet generation on the wall surface. The temperature along the liquid- 
vapor interface does not change significantly according to the Gibs-Thomson relation (see Shyy 
1994a), because the curvature of the interface is moderate The liquid-vapor interface can be 
reasonably assumed as isothermal, consequently the thermal capillary effect can be neglected 
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Hence the phase-change process is entirely driven by the internal heat generation and heat 
conduction. The rate of evaporation or condensation is determined by the local temperature 
distribution. The time dependent and non-uniform internal heat generation results in a transient 
temperature distribution and causes the motion and deformation of the interface. 

The internal energy formulation of the energy equation for this phase-change process with 
the constraint of constant volume is used The weak form of the energy equation for conduction is 

P% = V(kW) + 0,„ (2.9) 

at 

The source term O gen is the internal heat generation rate, representing the thermal power 
density released via nuclear fission reaction, which varies with time and space according to the 
following relation: 


Q gen = pG(0-sin(y) 


( 2 . 10 ) 


where G (t) is a time dependent volumetric heat generation rate, p is the local phase density, I is 
the length of the container. Some additional closure relations, the same as Eqs.(2.7) and (2.8) are 
required for the multiphase condition in the container. 


The governing equation (2.9) is normalized and expressed for different phases as: 

/ = 0 


p l ^=*yo+Q gtn , 

do _ , - _ 1 df 

p f~3r = KfV 6 + ® genJ ~ Pf St & x 


0</<l 
/ = ! 


( 211 ) 


where S, is the Stefan number defined as c v T s ,(/Ae. The nondimensional heat generation rate is 
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Qgen, = P,G(T)s\n(7lX) 


( 2 . 12 ) 


where G(z) =G(t)ftr/c v T s0 ). The initial saturation temperature T s ,o and the length of the container 
are chosen as the characteristic temperature scale and length scale, respectively. The characteristic 
time scale is tr=L 2 /a v which is the thermal diffusion time scale in the vapor phase. The E-based 
update method discussed previously is used for the iterative update of the vapor fraction. 


2.3 Thermophysical Properties of Uranium Tetrafluoride 

The nuclear fuel uranium tetrafluoride, UF 4 , is used to provide heat to the system. The 
thermophysical properties of its liquid and vapor at a reference state of T 0 =2000K and, P 0 =lMpa 
(10 atm ), are given in Table 2. 1 according to Watanabe and Anghaie (1993) and Anghaie (1992). 
Table 2. 1 Thermophysical Properties of UF 4 Liquid and Vapor at T 0 =2000K, P 0 =lMpa. 


Properties 

Cp 

Cv 

R 

P 

K 

P 

Pr 

Unit 

i/iWiil 

EMM 

EMM 

ES9 


BMW Ml 


Vapor 

345 

320 

25 

20 

0.06 

Effli 

0.53 


480 

480 


6000 

0.5 

naifgfil 

0.89 


where R is the gas constant. Some assumptions and simplifications have to be made based on the 
order and scaling analysis as discussed before. All the properties, except the vapor density, are 
assumed constant. The vapor phase is assumed to be the an ideal gas The vapor density is 
obtained from the equation of state for a perfect gas. In the governing equation (2.11), the 

— K i 

thermal conductivity ratio of liquid to vapor is k - — = 8.3 , and the density ratios of liquid to 

K v 

_ p, 

vapor at the referential state isp 0 = = 3x10 , where the subscript 0 denotes a referential 

A\o 

state. The Stefan number, St=CvT 0 /Ae, is 0.67 where Ae is the latent heat which is involved with 
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phase change. The characteristic time scale tr=L 2 /a v =280 seconds. The dimensionless time 
dependent parameter for heat generation rate expressed in Eq.(2. 12) is given by 


G(r) = 4r 0<r<1.5 

G(r) = 4 x15=6 >1.5 


(2.13) 


Thus the real parameter B(t) is 


G(t) = 33 xt 

G(0 = 33x420 =13.86 


(kW/kg) 


0 < t < 420 
t > 420 


(2.14) 


Correspondingly, the maximum thermal power density in the liquid phase, at the final 
steady state or the full power state when t > 420 s, is Q gen l | max = p,G(t > 4205) = 83 (MW ! m 3 ). 
The average power density, including the spatial distribution and vapor phase, for the entire fuel 
cell is 

Q g en. Lr age = 05 x Q genl | max ft sm(ar>fc = 0.5 X Q gen I | max x^ = 26.5 (MW / m 3 ) (2. 1 5) 


2.4 Results and Discussion 

The second order central difference scheme for the spatial differential terms and the fully 
implicit marching scheme for temporal term are adopted to solve Eq. (2. 1 1). A uniform mesh of 
101x41 grids covers the liquid phase and vapor phase within the container, and another non 
uniform mesh with 101x5 grids covers the side wall region. The outer surface of the side wall is 
cooled by a constant heat flux. 

Because of symmetry, only the right half of the container is considered Initially, the lower 
half of the container is filled with liquid and the upper half is occupied by vapor The aspect ratio 
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of the container is 2.5. The initial saturated vapor and liquid temperature is 2000K. At t=0, the 
side wall is cooled at a constant heat flux, and the internal heat starts to be generated. 

The evolution of the evaporation/condensation process is shown in Fig.2.1. Using the 
phase identifier which is the local vapor phase fraction, the distribution of the liquid and vapor 
phases and the position of the interface are determined. The topology of the liquid-vapor interface 
and the temperature distribution at various time instants are displayed in the top and bottom rows 
of Fig.2.1, respectively. The internal heat generation evaporates the liquid, while the side wall 
cooling tends to condense the vapor. Globally, the evaporation of liquid and condensation of 
vapor take place simultaneously, as shown schematically in Fig. 1.2. The local phase-change 
process, whether evaporation or condensation, is dominated by the local internal heating or by the 
wall cooling effect, since the internal heat generation is a spatial function. The evaporation is 
dominant near the central axial position where the heat generation intensity is higher and 
evaporates more liquid. The condensation caused by the side wall cooling is dominant near the 
axial end and side wall where the heat generation density is lower. The simultaneous evaporation 
and condensation cause the motion and deformation of the interface. Because of the time 
dependent heat generation, the local evaporation and condensation, the motion and deformation 
of the interface continue. 

As the process evolves, certain amount of liquid is evaporated into vapor due to the 
internal heat generation which is increasing at the early period, as Eq.(2.14) indicates. Meanwhile, 
some amount of vapor near the side wall is condensed into liquid due to the cooling of the side 
wall. During earlier periods, a short piece of liquid film is formed along the side wall owing to the 
condensation, and the liquid-vapor interface becomes concave because of the evaporation. With 
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time increasing, the interface keeps moving and deforming. The liquid film grows longer and 
thicker, and the position of the interface at the centerline of the container descends lower and 
lower. When t = 1 12 s, or r = 0.4, the liquid film grows up to the top end and covers the entire 
side wall. The liquid phase which is initially filled in the lower part of the container now forms a 
long film covering and residing along the entire side wall. At t = 392 s ( r=1.4 ), the interface at 
the low central side descends quite close to the bottom end. After t = 420 s ( t=J.5 '), the internal 
heat generation rate is a constant as indicated by Eq. (2.14) and the cooling heat flux at the side 
wall is still kept constant, the system reaches a steady state and the interface no longer moves. 
The topology of the interface strongly depends on the temperature distribution and the internal 
heat generation rate distribution which now is independent of time but only a spatial function. The 
shape of the interface follows approximately the heat generation distribution in a form of chopped 
sine function. 

During the transient process, the temperature distribution varies with the time. The 
temperature contours for the entire system become denser and denser gradually, since the internal 
heat generation is increasing. After the steady state is reached, the axial and radial temperature 
distributions follow approximately the forms of chopped sine and cosine function, respectively, 
due to the nature of the heat generation distribution as shown in Eq. (2.12). There is a 
temperature gradient jump at the interface due to the release or absorption of the latent heat 
associated with the phase change. The discontinuity of the temperature gradient can be observed 
more clearly in a three dimensional temperature distribution displayed in Fig.2.2. There is an 
obvious fold on the temperature surf where the liquid-vapor interface is located exactly The 
maximum temperature near the central axial position in the vapor region is about 4000K. The 
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axial temperature profile on the side wall is also in an approximate form of chopped sine function, 
but much flatter than that in the vapor region. The liquid is a much better conductor than the 
vapor, the liquid film covering the side wall is more isothermalized. This isothermalization effect is 
expected to satisfy the design feature of the multiphase fuel cell to be possibly applied in space. 
The formation of the liquid film and its full coverage on the side wall can prevent the side wall 
from being overheated, which is also the expectation for the application of such multiphase fuel 
cell. 

The transient characteristic of this system can also be seen from the variation of the 
saturation temperature with the time shown in Fig. 2.3. The saturation temperature decreases a 
little bit in the early time period when the internal heating is very weak, the side wall cooling 
effect and condensation are dominant. As the internal heating rate increases, the evaporation 
becomes dominant and the saturation temperature rises. After t = 420 s ( r = 1.5 ), as the system 
gradually reaches the steady state, the evaporation and condensation are balanced and the 
saturation temperature reaches a constant about 3000K. 

The operation in the zero gravity environment is ideal and unrealistic. The study 
conducted in this chapter, however, provides a basis for further discovering and exploring the 
fundamental performances and features of such multiphase and internally heated system. The 
present results will be used as reference to be compared with the performances and features of 
this system operated under micro-gravity conditions, which will be presented in the next chapter. 
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Figure 2.2 Temperature surf of the system at steady state ( 0-g ) 

( r = RxL = 0-0.01 m, x=XxL = 0 - 0.05 m, the real temperature is 1=0x10=0 - 4000K) 



Figure 2.3 Variation of saturation temperature ( 0-g ) 

(the real time period: 0 - 2800 s , the real temperature : T=0xT o =0x2OOOK ) 
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Chapter m. Modeling of Thermal Performance of Multiphase UF 4 Fuel Element 

at Micro-Gravity and Normal Gravity 

3.1 Brief Preview 

So far, only conduction dominated effects under limited conditions have been taken into 
account for the bulk evaporation and condensation processes. Like all phase change processes, 
evaporation and condensation are also inevitably subject to the influences of convection induced 
by either buoyancy or surface tension variations The existence of convection can substantially 
affect the flow field, the temperature field, and subsequently the performance of the phase change 
process. 

Numerous investigations on the phase change with convection have been conducted, 
which can be found in Ostrach (1983), Crank (1984), Brice (1986), Brown (1988), Gau and 
Viskanta (1986), Nadarajah and Narayanan (1990), Shyy et al. (1990, 1991, 1992a, b, 1994a, 
1995) and Shyy (1994a, b, 1995), etc. 

In this chapter, the modeling of the bulk evaporation and condensation processes, 
associated with internal heat generation and under various gravity levels, is presented. The system 
to be modeled is Fig. 1.2. The nuclear fuel uranium tetrafluoride, UF4, is used to provide heat to 
the system. The top wall and bottom wall are insulated. The side wall is cooled externally by a 
constant outward heat flux which removes heat from the system and tends to condense the vapor 
phase. The internal heat generation provides heat to evaporate the liquid phase. Since the system 
is to be operated under different gravity conditions, there is convection, both in the liquid and 
vapor phases, induced by buoyancy force due to the density variation. The temperature along the 
liquid-vapor interface will not change significantly according the Gibs-Thomson relation (see Shyy 
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1994a, 1995), because the curvature of the interface is moderate. The liquid-vapor interface can 
be reasonably assumed as isothermal, thus the thermal capillary effect can be neglected. The phase 
change processes are controlled by conduction and convection. The performances of the phase 
change can be quite different at various gravity condition, as we can expect The flow field, the 
temperature distribution and the motion of the interface are strongly dependent on the g-levels, as 
will be demonstrated later. 

3.2 The Governing Equations 

The governing equations describing the phase change process with the natural convection 
in both phases, in the cylindrical coordinate system, in original dimensional form are 

% + V*(aF) = 0 (3 1) 

at 

^ + V.(p,K») = -^ + V.( ft Vu)- p,g (3.2) 

at at 

— + V»(p,i 7 v) = -^ + V«(//,Vv) (3.3) 

dt 

^^ + V.(p,Fc v 70=^V 2 7+ A G(Osin(^y)-Ae^^-A e V.( A ^) (3.4) 

at L at 

where the subscript /=/, v denotes the liquid and vapor phases respectively, p v ,o is the referential 

vapor density, G(t) is the volumetric internal heat generation rate. The thermal expansion 

coefficient, /?, is introduced into the u-momentum equation to expand the buoyancy term, 

according to the Boussinesq approximation 

A=A..[i-0<r-r o )] OS) 
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Therefore we have 


~ Pig 

= gp ifi [fii(T-T 0 )-\] (3.6) 

= gp u oPM 0 ~ l )-gp,.0 

where the first term represents just the conventional buoyancy effect caused by the density 
variation in a single phase, while the second one is the global body force, or the bulk gravity force 
acted on the fluid at referential state. 

Equations (3.1) ~ (3 4) can be normalized as 

^ + V.(p,F) = 0 (3.7) 

or 

+ V . ( p,VU) = - Jp + Pr, V . (ft VU) + gRa. Pr. (9 - 1) - -pf\ (3 8) 

^f + V.(p,FV)=-^ + Pr.V.^V) (3.9) 

or cK 

^+W.(pfe)=Kyd + p i G(r)sm(7iX)-j i ^^-^V*(pff) (3.10) 

where g is the gravity factor representing the gravity level. The Rayleigh number, based on the 
ground gravity level and vapor properties, is defined 

Ra v = (3.11) 

®v,0^v,0 

The Prandtl number for vapor phase is 

Pr v = — ^ (312) 

«„n 
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The nondimensional internal heat generation rate is 


G(t) 


Gjt) 
CJs , 0 / tr 


(3.13) 


where tr=L 2 /a v ,o is the thermal diffusion time scale. 

The liquid Rayleigh number and Prandtl number are also implicitly included in the 
governing equations through the multiplication of the vapor Rayleigh number to the 
nondimensional thermal expansion coefficient and density ratio, and that of vapor Prandtl number 


to the nondimensional viscosity ratio. 

The gravity factor in Eq.(3.8), g , represents the gravity level. Such expression in the 
equation leads to an explicit and clear representation for the gravity effect and convenient 
adjustment for the computation. 

In obtaining the above equations, the following dimensionless variables are used 




P = 


P 

Pv.o(«v.O / Lf 



t • t 


(3.14) 


U = 


u 

Ko TTy 


v = 


(«v,0 ! L ) 


y = Ue x +Ve r 


where u r =a v .(/L is the thermal diffusion velocity scale in the vapor phase, e x and e r are the axial 
and radial unit vectors, respectively. Other velocity scales, such as buoyancy velocity scale, etc., 
can also be used to yield different forms of the governing equations. 

The thermophysical properties of UF 4 ’s liquid phase and vapor phase, at a referential state 
of T 0 =2000K and, P 0 =lMpa (10 atm ), are given in the Table 2.1 according to Watanabe and 
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Anghaie (1993) and Anghaie (1992). The initial density ratio of liquid to vapor, p/p v ,o, is 300, and 
the initial conductivity ratio, k/k V i0 , is 8.3. The vapor Rayleigh number Ra v based on the vapor 
phase properties and the ground gravity level is 1.6xl0 6 , vapor Prandtl number Pr v is 0.53, and 
Stefan number St is 0.67. The same amount of thermal power as shown in Section 2.3 and Eq. 
(3. 12) is generated and provide internal heating to the system. 

3.3 Numerical Treatment 

The fully implicit time marching, with time step size of 10' 3 , is used for the temporal 
difference. The hybrid scheme is used for the spatial difference, where the first order upwind 
scheme for the convection terms and the second order central difference scheme are used 
accordingly. 

The computation is unstable and difficult to converge when only the central difference 
scheme is used. One reason is that there are the interactions between the natural convection cells 
in the liquid phase and vapor phase, and the interaction is particularly strong near the liquid-vapor 
interface, as will be shown later. These interactions make the computation unstable and easily 
divergent. Another reason is that there is a computational grid limitation of cell Reynolds number, 
ReAx, or cell Peclet number, PeAx When the Rayleigh number or Grashof number is considerably 
large, a very thin grid has to be adopted, which is time consuming and expensive. As we know, 
for example, the natural convection velocity scale is u r =(Gr v ) ,/2 v/L, Gr v is the Grashof number 
based on vapor properties. Thus the cell Reynolds number, Re^=u r Ax/L=(Gr v ) ,/2 AX. It is 
required that Re M be not greater than 2 for the computational stability. To satisfy this, if Gr v =10 6 , 
zLTmust be smaller than 1/500. Such a dense mesh is not realistic for the computation capacity of 
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current ordinary work stations. By adopting under-relaxation and more inner-loop iterations for 
convection equations (see Shyy 1994a), etc., higher Re&x, i.e, 10-15, can be tolerated for some 
simple problems of natural convection in a single phase system, such as those of cases in chapter 
5. Associated with the phase change and the interactions of convection cells between two phases, 
the computational stability becomes more difficult to be guaranteed. Very dense meshes, i.e., 
41x101, 81x201, were tested for the cases of g = 10~ 3 , Rav=10 6 , and the central difference 
scheme failed. 

Therefore, the first order upwind scheme have to be used for the convection terms to 
guarantee the computational stability and convergence. The price paid is the lower accuracy of the 
results. In the following computations, the mesh with 51x101 uniform grids is used. With such 
mesh and 2000 time steps, the typical CPU time in the Workstation of SunSparc 20 is about 10 ~ 
15 days. 

3.4 Results and Discussions 

The evaporation and condensation processes under micro-gravity and normal gravity 
conditions, g = 10" 3 , and 1, are modeled. The dependence of the phase change performance on 
the gravity level is shown clearly. In the following results, the evolution of the liquid-vapor 
interface topology, the temperature distribution, the flow pattern of natural convection in liquid 
and vapor phases, are demonstrated. 

At the initial state, the liquid and vapor are assumed to be subcooled and superheated, 
respectively. The dimensionless temperature ranges from 0.8 to 1.0 in the liquid phase and from 
1.0 to 1.2 in the vapor phase. The liquid and vapor are saturated at the interface with the initial 
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saturation temperature of 1.0. The aspect ratio is 2.5. Because of the. axial symmetry, only right 
half of the container with the liquid and vapor is modeled. The internal heat generation provides 
heat both in the liquid and vapor phases according to Eqs.(2.13), (2.14). At v=0, the side wall 
also starts to be cooled with an outward constant heat flux. The concurrent internal heating and 
external cooling lead to the simultaneous bulk evaporation and condensation. As time marching, 
more rind more heat is generated until t = 420 s (t= 1.5) when the heat generation keeps constant 

3.4.1 At the micro-gravity of 10' 3 -g 

The results of the evaporation/condensation process is shown in Figure 3 . 1 through Figure 
3.5. Using the phase identifier which is the local vapor phase fraction, the distribution of the liquid 
and vapor phases and the position of the interface are determined. The topology of the liquid- 
vapor interface and the temperature distribution at various time instants are displayed in Figures 

3.1 and 3.2, respectively. The flow fields in the vapor and liquid phases are presented in Figures 
3.5(a), (b), respectively. 

The internal heat generation evaporates the liquid, while the side wall cooling tends to 
condense the vapor. Globally, the evaporation of liquid and condensation of vapor take place 
simultaneously. The local phase-change process, whether evaporation or condensation, is 
dominated by either the local internal heating or the wall cooling effect, since the internal heat 
generation is a spatial function. The evaporation is dominant near the central axial position where 
the heat generation intensity is higher and evaporates more liquid. The condensation caused by the 
side wall cooling is dominant near the side wall. The simultaneous evaporation and condensation 
cause the motion and deformation of the interface. Because of the nature of time dependent heat 
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generation, the local evaporation and condensation, the motion and deformation of interface 
continue 

As the process evolves, certain amount of liquid is evaporated into vapor due to the 
internal heat generation which increases at the early time period, as Eqs.(2. 13) and (2.14) indicate. 
Meanwhile, some amount of vapor near the side wall is condensed into liquid due to the cooling 
of the side wall. At earlier periods, a short piece of liquid film is formed along the side wall owing 
to the condensation, and the liquid-vapor interface becomes concave because of evaporation 
With time increasing, the interface keeps moving and deforming. The liquid film grows longer and 
thicker, and the position of the interface at the centerline of the container descends lower and 
lower. When t = 280 s or z =1.0, the liquid film grows up to the top end and covers the entire 
side wall. At t = 448 s or r=1.6, the interface at the low central side descends to the lowest point. 
The liquid phase, which initially fills in the lower part of the container, now forms a long film 
covering and residing along the entire side wall. After t = 420 s ( z= 1.5), the internal heat 
generation rate is a constant, as indicated by Eq. (2. 14), and the cooling heat flux at the side wall 
is also kept constant, the system reaches a steady state and the interface no longer moves The 
topology of the interface strongly depends on the temperature distribution and the internal heat 
generation rate distribution which now is independent of time and only a function of space. 

The shape of the interface does not follow the heat generation distribution in a form of 
chopped sine function, as it does in the case of zero gravity. The evolution of the interface 
topology is different from that of the case of zero gravity where no buoyancy convection is 
considered. Under the micro-gravity conditions, and due to the weak gravity effect, the buoyancy 
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induced convection is not strong but it does change the overall characteristics of the phase-change 


process. 

During the transient process, the temperature distribution varies with time The 
temperature contours for the entire system become denser and denser, since internal heat 
generation is increased. After steady state is reached, the axial temperature distribution does not 
follow exactly the form of chopped sine function, as it does in the case of zero gravity. There is a 
temperature gradient discontinuity at the interface due to the release or absorption of the latent 
heat associated with the phase change and the conductivity difference between vapor and liquid 
The discontinuity of the temperature gradient can be observed more clearly in the three 
dimensional temperature distribution displayed in Fig.3.3. There is an obvious fold on the 
temperature surf where the liquid-vapor interface is precisely located. The temperature gradient in 
the liquid region is greater than that in the vapor region, since the thermal power density is 
proportional to the mass density and more heat is generated in the liquid phase. The axial 
temperature profile on the side wall is much flatter than that in the vapor region Liquid is a much 
better conductor than vapor, the liquid film covering the side wall is more isothermalized. This 
isothermalization effect is expected to satisfy the design feature of the multiphase fuel cell for 
potential application in space. The formation of the liquid film and its full coverage on the side 
wall can prevent the side wall from being overheated, which is also the expectation for the 
application of such multiphase fuel cell The temperature distribution at each time instant is 
different from that in the case of zero gravity, since the convection induced by gravity, though 
very weak in this case, affects the temperature field. The convection are be discussed later. 
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The transient characteristic of this system can also be seen from the variation of the 
saturation temperature with the time shown in Fig. 3.4. The saturation temperature decreases a 
little bit in the early time period when the internal heating is very weak, the side wall cooling 
effect is dominant. As the internal heating rate increases, the evaporation becomes dominant and 
the saturation temperature rises. After t = 448 s (t=J.5), as the system gradually reaches the 
steady state, the evaporation and condensation are balanced and the saturation temperature 
maintains a constant. 

Figures 3.5(a) and (b) display the flow patterns of the natural convection at different time 
instants in the vapor phase and liquid phase, respectively. Vapor and liquid are two phases with 
largely different properties, and their flow and heat transfer characteristics are so different that 
need to be studied separately. In particular, due to the large mass density difference between two 
phases, the stream functions for vapor and liquid phases are calculated and plotted separately. 

During earlier times, several counter-rotating convection cells appear in the liquid phase. 
It is noticed that the convection in the liquid phase, owing to its higher density, is more intense 
than that in the vapor phase. The stronger convection cells in the liquid phase may induce 
convection in the vapor phase. On the other hand, the convection cells in the vapor phase also 
affect the flow in the liquid phase. These convection cells in both phases interact with each other 
and affect each other. Globally, the convection strength in both phases increases with time, since 
the temperature gradients across system become larger and larger gradually. 

It is found that the interaction between these convection cells also makes computation 
more difficult, and the strong under-relaxation and a lot of iterations, i.e., 500 ~ 1000, in each 
time step have to be used to guarantee the computational stability and convergence. 
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The global body force, -gp,, 0 in Eq.(3.6) or -gRa v Pr v p i 0 (——) in Eq.(3.8), always 

Pv * 0 

tends to keep the mass stay or move downward. The liquid phase - of larger mass density is 
subjected to greater body force than the vapor phase. The vapor phase will naturally occupy the 
space well above the liquid phase in case no bulk condensation involved. In this case, when the 
condensation occurs along the side wall due to its cooling, some amount of liquid is accumulated 
at the side wall. The continuing condensation and the accumulation of liquid lead to the formation 
of a thin and long liquid film along the side wall. Some liquid moves upward and downward, 
forming the natural convection due to the buoyancy effect caused by the density gradient and 
temperature gradient across the liquid film. Partial amounts of the liquid in the film is drawn by 
the body force downward and back to the liquid pool. A reflux cycle is formed by the phase 
transformation and the mass transfer. 

Compared to the results of the zero gravity case, the process of the liquid film formation 
and development is slower and less amount of liquid is accumulated in the film for the present 
case of micro-gravity. It is because of the gravitational force and the reflux effect. With the action 
of the gravity force, the liquid can not move upward as much and freely as it in the 0-g case. 
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Figure 3.3 Temperature surf and phase distribution at steady state ( 10' 3 -g ) 

( r = RxL = 0-0.01 m, x=XxL = 0 - 0.05 m, the real temperature is T=0xT o =O - 4000K) 



Time 


Figure 3.4 Variation of saturation temperature ( 10° -g ) 

(the real time period: 0 - 2800 s; the real temperature: T=0xT o =0x2OOOK ) 
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3.4.2 At the normal gravity of 1-g 

All the boundary and initial conditions are maintained the same as the case of 10 -g. The 
gravity factor is 1 in Eq.(3.8), The results are shown in Figure 3.6 through Figure 3.10. The 
topology of the liquid-vapor interface and the temperature distribution at various time instants are 
displayed in Figures 3.6 and 3.7, respectively. The flow fields in the vapor and liquid phases are 
presented in Figure 3.10. 

The global evolution of the phase change process is similar to the 10°-g case. As the time 
lapses, a piece of thin and long liquid film is formed and more and more liquid is accumulated 
along the side wall due to the cooling effect. At about t = 336 s ( x =1.2 ), the liquid film grows 
up to the top end and covers the entire surface of the side wall. After t =420 s (t=1.5), the system 
gradually reaches to a steady state when the shape of the liquid-vapor interface, the temperature 
distribution and the convection flow patterns do not change and evolve with the time. 

Significant differences appear in this normal gravity case, compared to the 10 3 -g case. 
First of all, the liquid film grows more slowly than it does with 10° -g. The liquid film is also 
thinner. Since the gravity effect is greater in the present case, the liquid film is more readily falling 
down. More liquid is drawn from the film back to the liquid pool and less liquid is accumulated in 
the film. And consequently the position of the horizontal surface of the liquid pool is higher than 
that in 1 0" 3 -g case. Secondary, the horizontal surface of the liquid pool becomes wavy, 
particularly in the position near the end of the liquid film. Such phenomenon does not appear in 
the case of 10' 3 -g. The liquid film, though thinner, carries more momentum owing to the greater 
gravity effect. It impacts the liquid pool and distorts the pool surface 
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Furthermore, the temperature distribution is quite different from those at 0-g and 10' 3 -g 
The temperature gradient across the liquid phase is much larger owing to the much intensive 
internal heat generation rate in the liquid phase. In particular, the temperature gradient across the 
liquid film is larger than that across the vapor phase, since the liquid film is thinner. Such a thin 
liquid film with large temperature gradient is a very good conductor for heat transfer. The 
internally generated heat in the vapor phase, which is less intensive than in the liquid phase, is 
barely accumulated and easier to be removed outward. Therefore the temperature gradient across 
the vapor region is nearly zero. 

Figure 3.10 shows the flow fields in the vapor and liquid phases. Compared to Figure 3.5, 
the flow patterns are different. The convection intensity is much larger in 1-g than that in 10' 3 -g, 
due to the larger gravity effect. The natural convection cells, in which some a mount of liquid 
moves upward and downward, are still formed in the thinner liquid film. The reflux of some 
amount of liquid back to the liquid pool can be seen more clearly. 
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t = 56 (s) 1 12 (s) 168 (s) 224 (s) 

x = 0.2 0.4 0.6 0.8 



t = 280 (s) 336 (s) 392 (s) 448 (s) 

x = 1.0 1.2 1.4 1.6 

Figure 3.7 The evolution of isotherm ( 1 -g ) 

( A0 = 0 02, AT=40K, the real temperature is T=0xT o =0x2OOOK ) 
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Figure 3.8 Temperature surf at steady state ( 1 -g ) 

( r = RxL = 0-0.01 m, x=XxL = 0 - 0.05 m, the real temperature is T=0xT o =O - 4000K) 



Figure 3.9 Variation of saturation temperature ( 10' 3 -g ) 
(the real time period. 0 - 2800 s, the real temperature: T=0xT o =0x2OOOK ) 
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The differences and the dependence of the phase change performance on the different 
gravity levels can be observed more clearly in Figure 3.11 where the topology of the liquid-vapor 
interface, the temperature distribution at the final steady state under three cases, 0-g, 10 -g and 1- 
g, are plotted together. The gravity influences on the thermal performance of the bulk liquid- 
vapor phase change process in the multiphase nuclear fuel cell are displayed explicitly Under all 
three gravity conditions, the liquid film is formed and covers the entire side wall, and the liquid 
film covering the side wall is more isothermalized at the wall surface, which can prevent the side 
wall from being overheated. At 0-g, the performance is controlled only by the heat transfer, such 
as internal heat generation and side wall cooling. At non-zero gravity conditions, its performance 
is controlled by both heat transfer and gravity effects. As gravity increases, the liquid film is 
thinner, the temperature gradient is larger across the liquid film and is smaller across the vapor 

phase. 

The model developed in this work can be used to provide very useful, and at least 
qualitative, information for the design and application of multi-phase nuclear fuel elements to be 
operated under different gravity conditions, particularly when the experiments at micro-gravity 
environment are expensive and not easy to be conducted. 

The current model can also be applied to some two-phase flow problems with or without 
phase change, for example, the film condensation, the liquid-vapor or liquid-gas flow problems in 
heat exchangers and material processing equipment. 
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(a) The topology of the liquid-vapor interface at the steady state 



Figure 3.11 Comparison of the interface topology and temperature distribution 

at different gravity levels 







Chapter IV. Conclusions 


A unique numerical method has been developed to model the dynamic processes of bulk 
evaporation and condensation processes, associated with internal heat generation and natural 
convection under different gravity levels The internal energy formulation, for the bulk liquid- 
vapor phase change problems in an encapsulated container, was employed. The equations, 
governing the conservation of mass, momentum and energy for both phases involved in phase 
change, were solved. 

The thermal performance of a multiphase uranium tetrafluoride fuel element under zero 
gravity, micro-gravity and normal gravity conditions was investigated. The evolution of the bulk 
liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and 
development of the liquid film covering the side wall surface, the temperature distribution and the 
convection flow field in the fuel element were obtained. Numerical modeling results show the 
strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity 
condition. The differences in the phase change performance of the multiphase fuel element at 
different gravity levels are quite significant. Under all three gravity conditions, 0-g, 10‘ 3 -g and 1-g, 
the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is 
more isothermalized at the wall surface, which can prevent the side wall from being overheated 
At 0-g, the thermal performance is controlled only by the heat transfer, such as internal heat 
generation and side wall cooling At non-zero gravity conditions, its thermal performance is 
controlled by gravity effect, as well as by heat transfer. As the gravity increases, the liquid film is 
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thinner, the temperature gradient is larger across the liquid film and smaller across the vapor 
phase. 

This investigation provides valuable, and at least qualitative, information for potential 
design and application of multi-phase nuclear fuel elements to be operated under different gravity 
conditions, particularly when the experiments at micro-gravity environment are expensive and not 
easy to be conducted. 

Further effort is worthful and recommended to extend the current model for liquid-vapor 
phase change problems to a unified triple-phase model. All phase change processes between these 
three phases, i.e., solid-liquid, liquid-vapor, solid- vapor, or solid-liquid-vapor, can be simulated 
This will contribute a lot for the research on materials processing. The current model can also be 
applied to some two-phase flow problems with or without phase change, for example, the film 
condensation, the liquid-vapor or liquid-gas flow problems in heat exchangers and material 
processing equipments. 
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